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Hydrogeology  and  Analysis  of  Aquifer  Characteristics  in 
West-Central  Pinellas  County,  Florida 

By  James  C.  Broska  ancf  Holly  L.  Barnette 


ABSTRACT 

The  U.S.  Geological  Survey,  in  cooperation 
with  Pinellas  County,  Florida,  conducted  an  inves¬ 
tigation  to  describe  the  hydrogeology  and  analyze 
the  aquifer  characteristics  in  west-central  Pinellas 
County.  A  production  test  well  and  four  monitor 
wells  were  constructed  in  Pinellas  County  at  Wals- 
ingham  Park  during  1996-97.  Water-quality  sam¬ 
pling,  static  and  dynamic  borehole  geophysical 
surveys,  and  hydraulic  tests  were  conducted  at  the 
wells  to  delineate  the  hydrogeology  at  Walsing- 
ham  Park.  A  9-day  aquifer  test  was  conducted  to 
determine  the  hydraulic  characteristics  of  the  aqui¬ 
fer  system  and  observe  the  changes  in  water  qual¬ 
ity  due  to  pumping. 

A  numerical  model  was  constructed  to  simu¬ 
late  the  aquifer  test  and  calculate  values  for 
hydraulic  conductivity  and  storage  coefficient  for 
permeable  zones  and  confining  units  at  Walsing- 
ham  Park.  Final  calibrated  values  for  hydraulic 
conductivity  for  the  different  permeable  zones  and 
confining  units  at  the  test  site  were  18  feet  per  day 
for  Upper  Zone  A,  750  feet  per  day  for  Lower 
Zone  A,  1  foot  per  day  for  Zone  B,  1x10'^  feet  per 
day  for  the  intermediate  confining  unit,  and  10  feet 
per  day  for  the  semiconfining  unit  separating 
Upper  Zone  A  and  Lower  Zone  A.  Final  calibrated 
values  for  storage  coefficient  were  S.lxlO'"^  for 
Upper  Zone  A,  8.6x10'^  for  Lower  Zone  A, 
2.6x10'^  for  Zone  B,  3. 1x10’"^  for  the  intermediate 
confining  unit,  and  4.3x10'^  for  the  semiconfining 
unit  separating  Upper  Zone  A  and  Lower  Zone  A. 
Estimates  of  transmissivity  for  Upper  Zone  A  and 
Lower  Zone  A  were  about  2,500  and  37,500  feet 
squared  per  day,  respectively. 


INTRODUCTION 

Pinellas  County,  located  in  west-central  Florida, 
has  limited  potable  water  supplies  and  is  dependent  on 
water  delivered  from  Tampa  Bay  Water  (formerly 
called  the  West  Coast  Regional  Water  Supply 
Authority)  which  withdraws  water  from  neighboring 
inland  counties  to  meet  tbe  water  demands  of  a  large 
population.  Although  the  population  of  Pinellas 
County  has  not  increased  significantly  since  1990, 
during  the  period  from  1970  to  1990  the  population  of 
Pinellas  County  increased  from  552,329  to  851,659 
(University  of  Florida,  1991).  In  1990,  the  county 
received  a  total  of  78  Mgal/d  of  water  from  Hillsbor¬ 
ough  and  Pasco  Counties  (Marella,  1992).  To  alleviate 
some  of  the  demand  for  this  delivered  water,  Pinellas 
County  is  exploring  the  possibility  of  using  the 
brackish- water  resources  within  the  Upper  Floridan 
aquifer  underlying  the  county. 

Brackish  waters  (water  with  dissolved-solids 
concentrations  between  1,000  and  10,000  mg/L)  are 
found  in  tbe  Upper  Floridan  aquifer  at  depths  ranging 
from  about  100  to  400  ft  below  land  surface  in  central 
Pinellas  County  (Cberry  and  others,  1970).  The  devel¬ 
opment  of  these  brackish-water  resources  through 
low-pressure  reverse  osmosis  could  provide  a  source 
of  potable  water.  Low-pressure  reverse  osmosis  is  a 
process  of  forcing  water  through  a  membrane  to 
remove  dissolved  solids,  thus  purifying  the  water; 
however,  the  process  is  economically  practical  only 
when  dissolved-solids  concentrations  are  less  than 
7,000  mg/L  (D.  Slonena,  Pinellas  County  Utilities, 
written  commun.,  1998). 

A  production  test  well  and  four  monitor  wells 
were  constructed  in  1996-97  at  Walsingham  Park  in 
west-central  Pinellas  County.  The  wells  were  sub¬ 
jected  to  various  geologic  investigation  methods  to 
determine  the  hydrogeology  at  the  test  site  and  the 
viability  of  brackish- water  development  in  the  area. 
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Purpose  and  Scope 

This  report  presents  the  results  of  a  study  hy  the 
U.S.  Geological  Survey  (USGS),  in  cooperation  with 
Pinellas  County,  to  evaluate  the  hrackish-water 
resources  underlying  the  Walsingham  Park  study  area 
in  west-central  Pinellas  County.  Results  from  the  field 
investigation  and  numerical  simulation  to  determine 
aquifer  characteristics  at  the  Walsingham  Park  test  site 
could  provide  valuable  information  for  future  develop¬ 
ment  of  hrackish  water  in  Pinellas  County  and  in  other 
coastal  counties  in  west-central  Florida.  The  report 
includes:  descriptions  of  the  hydrogeologic  framework 
of  the  study  site,  hydraulic  characteristics  of  the  aqui¬ 
fer  system,  and  description  of  the  ground-water 
quality. 

The  objectives  of  this  report  are  to 

1 . )  describe  the  hydrogeology  at  the  test  site,  and 

2. )  evaluate  an  aquifer  test  conducted  in  the  Upper 

Floridan  aquifer  at  the  test  site. 

Description  of  the  Study  Area 

Pinellas  County  is  a  peninsula  in  west-central 
Florida  bounded  by  Tampa  Bay  to  the  east  and  the 
Gulf  of  Mexico  to  the  west  (fig.  1).  The  county  is  part 
of  the  Gulf  Coastal  Lowlands  physiographic  region 
described  by  White  (1970)  which  consists  of  low 
angle  scarps  and  terraces  formed  during  several 
Pleistocene  sea-level  stands.  The  climate  of  Pinellas 
County  is  subtropical  with  an  average  rainfall  of 
53  inches  per  year  (in/yr).  Most  rainfall  occurs  during 
the  summer  months  between  June  and  September 
(Causseaux  and  Fretwell,  1983).  Land  surface  eleva¬ 
tion  at  the  Walsingham  Park  test  site  is  about  55  ft 
above  sea  level. 

HYDROGEOLOGIC  FRAMEWORK 

The  geology  of  Florida  consists  of  a  pre-Meso¬ 
zoic  igneous  and  metamorphic  basement  and  an  over- 
lying  thick  sequence  of  sedimentary  units  deposited 
during  transgressive  and  regressive  episodes  (Heath 
and  Smith,  1954).  The  hydrogeologic  framework 
underlying  the  study  area  consists  of  Eocene  or 
younger-age  rocks  and  makes  up  the  surficial  aquifer 
system,  the  intermediate  confining  unit,  and  the  Flori¬ 
dan  aquifer  system  (fig.  2).  This  multilayered  system 
consists  of  several  permeable  zones  interbedded  with 
zones  of  variable  permeability  that  function  as  either 
confining  units  or  semiconfining  units.  A  semiconfin- 


ing  unit  in  this  report  refers  to  lower-permeability 
carbonate  units  that  do  not  produce  large  quantities  of 
water  to  wells  when  pumped  (Knochenmus  and  Swen¬ 
son,  1996).  The  hydrogeologic  information  compiled 
from  wells  at  the  test  site  was  correlated  with  hydro¬ 
geologic  information  from  previous  studies  in  west- 
central  Pinellas  County.  A  generalized  hydrogeologic 
section,  oriented  north-south,  was  constructed  from 
natural-gamma  geophysical  logs  (fig.  3). 

Surficial  Aquifer  System 

The  surficial  aquifer  system  (SAS)  in  Pinellas 
County  is  the  shallowest  water-bearing  formation  in 
the  study  area  and  consists  of  undifferentiated  sands 
and  clays  which  change  in  composition  laterally  and 
vertically  (Causseaux,  1985).  The  post-Miocene  age 
surface  deposits  in  Pinellas  County  lie  unconformably 
above  the  Hawthorn  Group.  Surface  deposits  are  char¬ 
acterized  by  a  base  of  marginal  marine  beds  overlain 
by  sandy,  marine  terrace  deposits  and  capped  by  aeo- 
lian  sand  deposits.  The  thickness  of  the  SAS,  which 
ranges  from  a  few  feet  to  50  ft,  restricts  the  amount  of 
potable  water  the  aquifer  supplies;  however,  the  non- 
potable  water  is  utilized  for  irrigation  (Gilboy,  1985). 
The  SAS  at  the  Walsingham  Park  site  is  about  35  ft 
thick.  Throughout  Pinellas  County,  the  water  table  is 
less  than  15  ft  below  land  surface  during  dry  periods 
and  close  to  land  surface  during  wet  periods.  Precipi¬ 
tation  is  the  main  source  of  recharge. 

Hydraulic  properties  of  the  SAS  are  highly  vari¬ 
able  due  to  the  diversity  in  physical  characteristics  of 
the  sediments.  Grain  size,  sorting,  and  thickness  vary 
throughout  the  system  and  greatly  influence  the  per¬ 
meability  and  porosity  of  the  system.  Values  for  spe¬ 
cific  yield  range  from  0.1  to  0.35  (Causseaux,  1985). 
Horizontal  hydraulic  conductivity  is  estimated  to 
range  from  13  to  33  ft/d  (Cherry  and  Brown,  1974; 
Sinclair,  1974)  and  vertical  hydraulic  conductivity  is 
estimated  to  range  from  0.36  to  13  ft/d  (Sinclair,  1974; 
Hutchinson  and  Stewart,  1978). 

Intermediate  Confining  Unit 

The  surficial  and  Floridan  aquifer  systems  are 
separated  by  the  low-permeability  intermediate  con¬ 
fining  unit  (ICU).  The  ICU  occurs  within  the  upper 
Miocene  undifferentiated  Arcadia  Formation  of  the 
Hawthorn  Group.  The  formation  consists  of  interbedded 
quartz  sands,  clays,  and  carbonates.  The  clay  beds  in 
the  Arcadia  Formation  act  as  confining  beds 
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82°50' 


82°45' 


Albers  Equal-Area  Conic  projection 

Standard  Parallels  29°30'  and  45°30',  central  meridian  -83°00' 


Figure  1.  Location  of  the  study  area  in  west-central  Pinellas  County,  well  names,  and  hydrogeologlc 
section  A-JK. 


throughout  Pinellas  County  (Heath  and  Smith,  1954). 
The  thickness  of  the  ICU  averages  about  90  ft  and 
ranges  from  50  to  140  ft.  At  the  Walsingham  Park  site, 
the  ICU  is  about  100  ft  thick.  Due  to  the  heteroge¬ 
neous  nature  of  the  unit,  a  wide  range  of  hydraulic 
properties  has  been  reported.  Vertical  hydraulic  con¬ 


ductivity  for  the  ICU  ranges  from  l.SxlO  '^  to  6.9x10'^ 
ft/d  (Sinclair,  1974;  Black,  Crow  and  Eidness,  Inc., 
1978).  Ueakance  values  reported  in  previous  investi¬ 
gations  range  from  1x10"^  to  1.5x10"^  (ft/d)/ft  (Black, 
Crow  and  Eidness,  Inc.,  1978;  Seaburn  and  Robertson, 
Inc.,  1983;  Brown  and  Associates,  Inc.,  1986). 
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site,  Pinellas  County,  and  equivalent  model  units  (modified  from  Knochenmus  and  Swenson,  1996). 
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Floridan  Aquifer  System 

The  Floridan  aquifer  system  consists  of  lime¬ 
stones  and  dolomites  and  includes  the  Upper  Floridan 
aquifer,  the  middle  confining  unit,  and  the  Lower 
Floridan  aquifer.  The  Upper  Floridan  aquifer  in  Pinel¬ 
las  County  is  composed  of  the  Avon  Park  Formation, 
the  Ocala  Limestone,  the  Suwannee  Limestone,  and 
the  Tampa  Member  of  the  Arcadia  Formation  (fig.  2). 
The  Avon  Park  Formation  consists  of  marine  lime¬ 
stones  interhedded  with  dolostones.  The  limestones 
are  predominantly  well  indurated  packstones  which 
are  highly  fossiliferous.  Overlying  the  Avon  Park  For¬ 
mation  is  the  Ocala  Limestone.  The  Ocala  Limestone 
consists  of  moderately  indurated  grainstones  and 
packstones.  The  foraminiferal  limestones  in  the  unit 
commonly  contain  henthic  foraminifera,  hryozoan, 
echinoid,  and  mollusk  fossils.  The  Ocala  Limestone 
has  a  gradational  contact  with  the  overlying  Suwannee 
Limestone.  The  Suwannee  Limestone  is  predomi¬ 
nantly  foraminiferal  packstone  to  grainstone  with 
abundant  mollusk  and  microfaunal  remains  (Miller, 
1986).  Overlying  the  Suwannee  Limestone  is  the 
Tampa  Member  of  the  Arcadia  Formation,  which  is 
the  lower  Miocene  unit  of  the  Hawthorn  Group 
(Miller,  1986).  The  Tampa  Member  lies  unconform- 
ably  on  the  Suwannee  Limestone  and  consists  of  lime¬ 
stone  with  dolostone,  sands,  and  clays  present  in 
minor  amounts.  Mollusks,  foraminifera,  and  algae  are 
the  common  fossil  molds  present  in  the  unit.  The  unit 
varies  from  mudstones  to  packstones,  but  is  predomi¬ 
nately  wackestones. 

In  Pinellas  County,  the  Upper  Floridan  aquifer  is 
subdivided  into  four  permeable  zones  generally  sepa¬ 
rated  by  semiconfining  units  as  defined  by  Hickey 
(1982).  The  zones  are  labeled  alphabetically  with 
increasing  depth  from  A  to  D.  The  focus  of  this  report 
is  the  uppermost  part  of  the  Upper  Floridan  aquifer, 
specifically  the  Tampa  Member  of  the  Arcadia  Forma¬ 
tion  of  the  Hawthorn  Group  (Zone  A)  and  the  Suwan¬ 
nee  Limestone  (Zones  A  and  B)  because  these  are  the 
only  zones  that  contain  brackish  water  at  concentra¬ 
tions  suitable  for  low-pressure  reverse  osmosis.  Fur¬ 
ther  discussions  will  be  restricted  to  Zones  A  and  B. 

Zone  A 

Zone  A,  the  shallowest  and  freshest  permeable 
zone  in  the  Upper  Floridan  aquifer,  is  made  up  of  the 
Tampa  Member  and  the  upper  part  of  the  Suwannee 
Limestone.  In  Pinellas  County,  the  thickness  of  Zone  A 
averages  180  ft  and  ranges  from  approximately  115  to 


250  ft.  Reported  transmissivities  for  Zone  A  through¬ 
out  Pinellas  County  range  from  10,000  to  40,000  ft^/d 
(Hickey,  1982;  Seaburn  and  Robertson,  Inc.,  1983; 
Brown  and  Associates,  Inc.,  1986).  Values  for  stora- 
tivity  range  from  4x10“^  to  8x10“^  (Hickey,  1982). 
Hickey  (1982)  divided  Zone  A  into  Upper  Zone  A 
(UZA)  and  Lower  Zone  A  (LZA),  separated  by  a 
discontinuous  semiconfining  unit.  At  the  Walsingham 
Park  site,  a  25 -ft  thick  semiconfining  unit  consisting  of 
an  abundance  of  marine  clay  separates  the  140  ft  thick 
UZA  from  LZA,  which  is  about  50  ft  thick  (fig.  2). 

Zone  B 

In  Pinellas  County,  Zone  B  underlies  Zone  A. 
Thickness  of  Zone  B  averages  about  150  ft  and  ranges 
from  about  125  to  170  ft.  The  upper  part  of  Zone  B 
functions  as  a  semiconfining  unit  in  some  parts  of 
Pinellas  County  and  is  considered  a  nonproducing 
zone  in  these  areas  (Hickey,  1982;  Knochenmus  and 
Swenson,  1996).  Reported  hydraulic  properties  for 
this  sequence  are  1.3x10  to  2  ft/d  for  vertical  hydrau¬ 
lic  conductivity,  0. 1  ft/d  for  horizontal  hydraulic  con¬ 
ductivity,  and  7x10"^  (ft/d)/ft  for  leakance  (Hickey, 
1982;  Seaburn  and  Robertson,  Inc.,  1983;  and  Brown 
and  Associates,  Inc.,  1986).  The  lower  part  of  Zone  B 
is  more  permeable  and  is  comprised  of  the  lower  part 
of  the  Suwannee  Limestone.  Thickness  of  this  section 
of  Zone  B  averages  about  62  ft  and  ranges  from  50  to 
75  ft.  Transmissivity  estimates  by  Hickey  (1982) 
suggest  a  value  of  about  5,000  ft  /d.  At  the  study  site. 
Zone  B  was  encountered  between  350  and  547  ft 
below  land  surface.  However,  no  apparent  changes  in 
permeability  in  the  upper  and  lower  parts  of  Zone  B 
could  be  determined  during  drilling. 

Delineation  of  Hydrogeologic  Units 
Underlying  the  Study  Area 

The  delineation  of  hydrogeologic  units  and  their 
associated  hydraulic  properties  at  the  Walsingham 
Park  test  site  was  accomplished  through  the  collection 
and  interpretation  of  lithologic  and  stratigraphic  data, 
water-quality  data,  borehole  geophysical  data,  and 
hydraulic  data.  During  the  well  construction,  rock  core 
description,  discrete  water-quality  sampling,  borehole 
geophysical  logging,  and  specific  capacity  tests  were 
performed.  Upon  completion  of  the  well  construction, 
a  9-day  aquifer  test  was  performed.  In  total,  one  test 
production  well  3  and  four  monitor  wells  were 
installed  at  the  test  site  (fig.  4). 
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NOTE:  Not  to  scale 


Figure  4.  Walsingham  Park  test  site,  including  location  of  wells  drilled  during  the  field 
investigation,  and  location  of  off-site  monitor  well  (general  location  is  shown  in  fig.  1). 


A  production  test  well  and  four  monitor  wells 
were  drilled  by  using  a  dual- wall  reverse-circulation 
rotary  drilling  method.  With  this  drilling  method,  pres¬ 
surized  air  is  forced  down  between  the  inner  and  outer 
barrels  of  the  drilling  stem.  Rock  and/or  sediment  is 
cut  by  the  drill  bit,  and  core  samples  rise  to  land  sur¬ 
face  through  the  inner  barrel.  The  first  well  drilled  was 
originally  designed  to  function  as  a  deep  exploratory 
well  and  was  drilled  to  a  depth  of  547  ft.  This  well  was 
drilled  to  define  the  geology,  water  quality,  and  verti¬ 
cal  distribution  of  producing  zones  and  confining  units 
at  the  test  site.  Data  collected  from  the  deep  explor¬ 
atory  well  were  used  to  determine  the  total  depths  and 
casing  depths  of  the  test  production  well  and  monitor 
wells.  After  drilling  and  data  collection,  the  deep 
exploratory  well  was  backplugged  to  a  depth  of  266  ft 
and  converted  into  the  Walsingham  Park  UZA  monitor 
well  1  (WPMW-1).  The  WPMW-1  and  all  other  moni¬ 
tor  wells  were  constructed  to  monitor  a  specific  pro¬ 
ducing  zone  or  lower  permeability  unit  to  evaluate  the 
responses  of  various  zones  to  pumping.  The  other 
wells  drilled  as  part  of  this  study  are  hereafter  referred 
to  as  the  test  production  well  1  (WPTPW-1),  the  LZA 
brackish-water  monitor  well  2  (WPMW-2),  the  inter¬ 
mediate  confining  unit  monitor  well  3  (WPMW-3), 


and  the  surficial  aquifer  system  monitor  well  4 
(WPMW-4).  Total  depths  and  casing  depths  of  all  the 
wells  constructed  during  this  study  are  provided  in 
figure  5. 

Water-Quality  Sampling 

Water-quality  samples  were  collected  by  the 
USGS  and  Pinellas  County  personnel  during  drilling  of 
WPMW- 1  to  determine  vertical  distribution  of  water 
quality  at  the  site.  Water  samples  were  collected  at  dis¬ 
crete  intervals  every  10  ft  during  drilling  and  tested 
onsite  for  specific  conductance,  temperature,  and  pH. 
When  abrupt  changes  in  field  measurements  occurred, 
water  samples  were  collected  and  analyzed  for  major 
ions  and  physical  properties  at  the  Pinellas  County 
Utilities  laboratory.  A  vertical  water-quality  profile 
showing  variations  in  dissolved  solids  with  depth  dur¬ 
ing  drilling  of  WPMW-1  is  presented  in  figure  6.  As 
illustrated,  the  dissolved-solids  concentration  begins 
increasing  approximately  350  ft  below  land  surface. 
Although  some  water  samples  collected  between  360 
and  425  ft  below  land  surface  contained  a  lower  con¬ 
centration  of  dissolved  solids,  this  anomaly  may  have 
been  caused  by  drilling  constraints  as  water  samples 
were  collected  over  a  period  of  3  consecutive  days. 
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LOCATION  OF  UNIT  BOUNDARIES  AT  DEPTH  BELOW  LAND  SURFACE,  IN  FEET 


EXPLANATION 


OPEN  BOREHOLE 


SCREENED  INTERVAL 


GROUT  OR  CEMENT 
WELL  CASING 


NOTE:  Before  backplugging  to  266  feet,  WPMW-1  was  referred  to  as  the  deep  exploratory  well. 


Figure  5.  Generalized  section  view  of  casing  depths  and  finai  depths  of  weiis  at  the  Waisingham 
Park  test  site. 
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DISSOLVED  SOLIDS, 

IN  MILLIGRAMS  PER  LITER 


Figure  6.  Relation  of  dissolved  solids  and  depth  from 
water  samples  collected  during  drilling  of  Walsingham  Park 
deep  exploratory  well  (WPMW-1). 


The  freshening  of  water  samples  between  drilling  days 
may  have  been  the  result  of  mixing  of  waters  in  the 
borehole.  However,  previous  investigations  in  northern 
Pinellas  County  have  shown  that  layering  of  water  with 
higher  dissolved-solids  concentrations  can  occur  over 
water  with  lower  dissolved-solids  concentrations 
(Knochenmus  and  Swenson,  1996).  The  nature  of  this 
water-quality  anomaly  at  the  Walsingham  Park  test  site 
could  not  be  determined  upon  completion  of  drilling. 

Additionally,  selected  samples  were  used  to 
show  the  correlation  between  specific  conductance  and 
dissolved-solids  concentrations  for  ground  water  repre¬ 
sentative  of  the  study  site  (fig.  7).  The  linear  regression 
closely  fits  a  straight  line  with  a  correlation  coefficient 
(r)  of  0.997.  The  correlation  of  the  two  parameters  was 
later  used  to  estimate  changes  in  dissolved-solids 
concentration  under  pumping  conditions;  however,  the 
relation  of  dissolved  solids  to  specific  conductance 
would  not  be  constant  for  the  entire  range  of  samples 
(Hem,  1985).  As  figure  7  illustrates,  the  slope  of  the 
linear  regression  is  different  for  specific  ranges  of 
conductance  values.  Changes  in  water  quality  that 
occurred  during  hydraulic  testing  are  discussed  further 
in  the  aquifer  test  section  of  this  report. 

Geophysical  Logging 

Borehole  geophysical  logging  was  used  to 
delineate  hydrogeologic  units,  flow  zones,  and  water- 
quality  changes  at  the  test  site.  Geophysical  logs  were 
collected  by  a  Pinellas  County  contractor  under  static 
and  pumping  conditions  on  WPMW-1.  A  lithologic 
log  from  WPMW- 1  was  also  compiled  during  drilling. 
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Figure  7.  Relation  between  dissolved  solids  and  specific 
conductance  of  selected  water  samples  collected  during 
drilling  of  Walsingham  Park  deep  exploratory  well  (WPMW-1). 


Logs  collected  during  static  conditions  included  cali¬ 
per,  natural-gamma  radiation,  fluid  resistivity,  and 
fluid  temperature.  Logs  collected  under  pumping  con¬ 
ditions  included  ascending  and  descending  flow  meter 
with  the  pump  placed  120  ft  below  land  surface. 

A  schematic  of  the  geophysical  logs  collected  from 
WPMW-1  is  shown  in  figure  8. 

Interpretations  of  the  geophysical  logs  in  con¬ 
junction  with  the  lithologic  log  were  used  to  determine 
the  distribution  of  flow  zones  and  semiconfining  units. 
Natural-gamma  radiation  logs  are  used  to  delineate 
clay  and  phosphate-bearing  deposits  by  recording  the 
naturally  occurring  radiation  coming  from  the  forma¬ 
tion  adjacent  to  the  borehole.  The  natural-gamma  log, 
when  correlated  with  the  WPMW-1  lithologic  log, 
confirmed  the  presence  of  clay  and  phosphate  deposits 
from  35  to  135  ft  below  land  surface.  This  sequence 
corresponded  to  the  ICU  at  the  study  site.  The 
decrease  in  the  natural-gamma  values  at  about  340  ft 
below  land  surface  indicates  a  region  of  low  phosphate 
content  within  the  Suwannee  Limestone. 

Flow  logs  measured  under  pumping  conditions 
are  used  in  conjunction  with  caliper  logs  to  delineate 
flow  zones.  Within  a  relatively  uniform  borehole 
diameter,  an  increase  in  counts  per  second  indicates  an 
increase  in  fluid  velocity,  thus  an  increase  in  flow. 

The  flow  meter  logs  collected  at  WPMW-1  indicate 
the  presence  of  a  major  flow  zone  from  around  300  to 
350  ft  below  land  surface,  corresponding  to  LZA. 


Hydrogeologic  Framework  9 


CALIPER 

GAMMA  RAY 

FLOW 

FLOW 

STATIC 

STATIC 

METER 

METER 

FLUID 

FLUID 

ASCENDING* 

DESCENDING* 

RESISTIVITY 

TEMPERATURE 

SURFACE 


74  77  80 


600 

1  INCHES  COUNTS  PER  SECOND  COUNTS  PER  SECOND  COUNTS  PER  SECOND 

OHM-METER 

DEGREES  FAHRENHEIT 

Figure  8. 

*  FLOW  METER  LOGS  CONDUCTED  AT  A  PUMPING  RATE  OF  387  GALLONS  PER  MINUTE 

Geophysical  logs  of  Walsingham  Park  deep  exploratory  well  (WPMW-1). 

Fluid  resistivity  and  temperature  logs  are  used 
to  identify  zones  of  water-quality  change.  The  fluid 
resistivity  and  fluid  temperature  logs  indicate  a  rapid 
water-quality  change  approximately  350  ft  helow  land 
surface  as  fluid  resistivity  decreases  and  fluid  tempera¬ 
ture  increases.  The  change  in  water  quality  at  350  ft 
helow  land  surface  was  noted  in  the  water-quality  pro¬ 
file  developed  from  water  samples  taken  during  drill¬ 
ing.  However,  unlike  the  water-quality  profile,  neither 
the  fluid  resistivity  nor  temperature  logs  suggest  a 
freshening  of  water  from  360  ft  helow  land  surface  to 
the  lowest  extent  of  the  logs  at  about  415  ft  helow  land 
surface  (fig.  8). 

Specific-Capacity  Tests 


The  values  were  estimated  hy  applying  specific- 
capacity  and  packer-test  data  to  empirical  equations 
developed  from  the  modified  Jacob  nonequilibrium 
equation  for  a  single-well  test  (Driscoll,  1986,  p.  1021). 
Due  to  the  difficulty  in  placing  packers  in  the  borehole 
and  leakage  problems  detected  around  the  packers, 
confidence  in  the  packer-test  data  was  low  for  deter¬ 
mining  hydraulic  parameters.  However,  UZA  from 
about  135  to  270  ft  below  land  surface  and  LZA  from 
about  300  to  350  ft  below  land  surface  were  deter¬ 
mined  to  be  the  primary  flow  zones  contributing 
ground  water  to  the  wells.  Although  LZA  had  a  higher 
specific  capacity,  it  was  located  directly  above  a  highly 
mineralized  water-quality  zone.  Therefore,  UZA  was 
more  suitable  for  a  longer-term  aquifer  test. 


Specific-capacity  tests  were  performed  on 
WPMW- 1  to  evaluate  the  various  hydrogeologic  units 
at  the  study  site.  Water-level  data,  collected  under  a 
range  of  pumping  rates,  were  analyzed  to  estimate 
hydraulic  parameters  of  the  system.  Specific-capacity 
tests  measure  well  yield  per  unit  of  drawdown  for 
specific  intervals  of  individual  hydrogeologic  units. 
Results  were  used  to  estimate  hydraulic  properties  of 
specific  intervals,  locate  permeable  zones,  and  deter¬ 
mine  the  well  and  casing  depths  of  the  monitor  and 
production  wells.  In  addition  to  specific-capacity  tests, 
packer  tests  were  performed  upon  completion  of 
WPMW-1  (prior  to  backplugging)  to  further  identify 
permeable  zones  at  the  site.  Transmissivity  and 
hydraulic  conductivity  values,  estimated  from  spe¬ 
cific-capacity  and  packer  tests  performed  on  the  deep 
exploratory  well  (WPMW-1),  are  listed  in  table  1. 


Aquifer  Test 

A  9-day  aquifer  test  was  conducted  at  the  test 
site  from  July  29  to  August  7,  1997,  to  estimate  the 
hydraulic  parameters  of  Zone  A  and  determine 
changes  in  water  quality  due  to  long-term  pumping. 
The  four  onsite  and  one  offsite  monitor  wells  were 
used  to  observe  responses  in  various  zones  to  pumping 
from  UZA  at  WPTPW-1.  Background  water-level  data 
collected  before  and  after  pumping  showed  that  head 
differences  between  zones  at  the  test  site  reflected  a 
downward  head  gradient  (fig.  9).  At  the  onset  of  the 
aquifer  test,  the  water  level  in  the  SAS  monitor 
(WPMW-4)  was  about  2.5  ft  above  the  water  level  in 
the  ICU  monitor  (WPMW-3),  whereas  the  water  level 
in  the  ICU  monitor  was  more  than  36  ft  higher  than 
the  water  level  in  the  UZA  monitor  (WPMW-1). 

Moreover,  the  water  level  in  the  UZA 


Table  1.  Results  from  specific  capacity  and  packer  tests  from  the  deep 
exploratory  well 

[ft,  feet;  fP/d,  cubic  feet  per  day;  (gal/min)/ft,  gallons  per  minute  per  foot;  fF/d, 
feet  squared  per  day;  ft/d,  feet  per  day] 


Depth  interval 
below  land 
surface 
(ft) 

Pump 

rate 

(fF/d) 

Total 

change 

In  head 
(ft) 

Specific 

capacity 

{gal/min)/ft 

Transmis¬ 

sivity 

(ft^/d) 

Hydraulic 

conductivity 

(ft/d) 

135-185 

11,744 

9.73 

6.3 

1,675 

34 

135-235 

11,934 

8.87 

7.0 

1,870 

19 

135-235 

59,864 

59.64 

5.2 

1,400 

14 

205-235 

9,336 

9.59 

5.1 

1,350 

45 

230-282“ 

17,901 

13.74 

6.8 

1,810 

35 

268-285“ 

9,182 

4.63 

10.3 

2,750 

162 

288-340“ 

18,289 

14.50 

6.6 

1,750 

34 

338-355“ 

9,624 

4.31 

11.6 

3,100 

182 

385-542“ 

9,124 

5.77 

8.2 

2,195 

33 

“Packer  test  interval. 


monitor  was  0.15  ft  higher  than  the  water 
level  in  the  LZA  monitor  (WPMW-2). 
The  test  production  well  (WPTPW-1)  was 
pumped  at  a  rate  of  220  gal/min  for  the 
length  of  the  test.  After  the  initiation  of 
pumping,  responses  were  noted  for  the 
various  monitor  wells.  WPMW-1 
responded  immediately  with  a  lowering 
of  water  level  in  the  first  minute  of  pump¬ 
ing.  WPMW-2  also  responded  quickly 
with  the  first  drawdown  noted  after  4 
minutes  of  pumping.  Contrasting  the 
changes  in  water  levels  at  these  two  wells, 
the  water  level  in  WPMW-1  decreased  by 
1.0  ft,  whereas  the  water  level  in  WPMW-2 
decreased  by  0.02  ft  after  4  minutes  of 
pumping.  Furthermore,  a  “reverse  water- 
level  fluctuation”  (Andreasen  and 
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WATER  LEVEL,  IN  FEET  ABOVE  SEA  LEVEL 


Figure  9.  Water  levels  In  all  monitor  wells  before,  during,  and  after  the  9-day  aquifer  test. 
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Brookhart,  1963)  was  observed  in  the  ICU  monitor 
well  WPMW-3.  The  water  level  in  the  ICU  began  to 
rise  immediately  after  pumping  began  and  declined 
when  pumping  ended  (fig.  9).  A  possible  explanation 
for  this  occurrence  is  that  the  release  of  water  from 
compressive  storage  produced  by  pumping  in  UZA 
caused  a  reduction  in  aquifer  volume  of  UZA,  which  in 
turn  caused  both  horizontal  and  vertical  deformations 
in  the  overlying  ICU.  These  deformations  of  the  ICU 
subsequently  induced  the  measured  water-level  rise. 
This  phenomenon  has  been  documented  and  is  often 
referred  to  as  the  Noordbergum  Effect  (Hsieh,  1996). 
No  response  due  to  pumping  was  denoted  in  the  surfi- 
cial  aquifer  system  monitor  well  WPMW-4.  Although 
water  levels  were  monitored  in  the  five  wells  for  the 
duration  of  the  test,  drawdowns  reached  steady-state 
within  the  first  3  hours  of  pumping  when  tidal  effects 
began  to  interfere  with  drawdown  measurements 
(fig-  9). 

Water-quality  data  were  collected  for  the  dura¬ 
tion  of  the  9-day  aquifer  test.  Specific  conductance 
probes  were  placed  in  WPMW- 1  and  WPMW-2  at 
depths  of  250  ft  and  330  ft,  respectively.  Discharge 
water  from  WPTPW-1  was  monitored  hourly  for  spe¬ 
cific  conductance  for  the  duration  of  the  aquifer  test. 
The  changes  in  specific  conductance  measured  during 
the  test  are  shown  in  figure  10.  Interestingly,  water- 
quality  changes  noted  in  the  UZA  and  the  UZA 
monitor  wells  exhibited  opposite  trends.  By  using  the 
relation  between  dissolved  solids  and  specific  conduc¬ 
tance  developed  earlier  (fig.  7),  an  inference  can  be 
made  regarding  the  changes  in  dissolved  solids  con¬ 
centrations  during  the  aquifer  test.  The  specific  con¬ 
ductance  in  WPMW-1  in  UZA  increased  from  about 
1,050  to  1,600  |lS/cm  corresponding  to  an  increase  in 
dissolved-  solids  concentration  of  about  372  mg/L. 

By  contrast,  the  specific  conductance  in  WPMW-2  in 
UZA  decreased  from  about  11,900  to  11,200  |lS/cm 
corresponding  to  a  decrease  in  dissolved-solids  con¬ 
centrations  of  about  440  mg/L  over  the  9  days  of  the 
test.  The  specific  conductance  of  the  composite  sam¬ 
ple  collected  from  WPTPW-1  increased  from  470  to 
488  p-S/cm.  Although  this  slight  increase  in  measured 
specific  conductance  can  be  considered  to  be  within 
the  error  of  measurement  of  the  conductance  probe, 
the  data  showed  a  consistent  trend  (fig.  10). 

The  opposite  trends  in  water-quality  changes 
noted  between  the  UZA  and  the  LZA  monitor  wells 
during  the  9-day  aquifer  test  are  difficult  to  explain. 


The  upward  hydraulic  gradient  between  LZA  and 
UZA  induced  during  the  aquifer  test  may  have 
resulted  in  upward  movement  of  more  saline  water. 
However,  the  freshening  of  water  in  the  LZA  monitor 
well  (WPMW-2)  during  the  aquifer  test  is  anomalous. 
As  figure  10  illustrates,  the  rate  of  change  of  measured 
specific  conductance  appears  to  decrease  from  the 
time  pumping  began.  It  is  possible  that  mineralized 
water  may  have  invaded  LZA  during  drilling  of  the 
deep  exploratory  well  and  may  have  evacuated  LZA 
during  the  aquifer  test,  thus  returning  to  native  condi¬ 
tions.  Since  the  aquifer  test  lasted  only  9  days,  no  data 
are  available  to  determine  the  specific  conductance 
values  at  longer  pumping  times. 

Data  and  the  results  of  the  aquifer  test  are  on  file 
and  available  at  the  USGS  Tampa  office.  Because  ana¬ 
lytical  methods  are  constrained  by  limitations  in  the 
simplification  of  a  flow  system  into  one-  or  two- 
dimensional  flow,  the  9-day  aquifer  test  was  analyzed 
by  using  numerical  methods.  Numerical  methods  can 
simultaneously  calculate  changes  in  head  and  advec- 
tive  transport  in  a  variable-density  flow  system  for 
multiple  observation  points. 

ANALYSIS  OF  AQUIFER 
CHARACTERISTICS 

A  numerical  modeling  approach  was  used  to 
simulate  the  effects  of  pumping  the  variable-density 
ground- water  flow  system  at  the  test  site.  The  model 
was  constructed  to  simulate  the  9-day  aquifer  test  con¬ 
ducted  at  the  site  and  to  refine  estimates  of  hydraulic 
parameters  such  as  hydraulic  conductivity  and  storage 
coefficient  determined  from  specific-capacity  tests  and 
previous  investigations. 

Numerical  Methods 

The  USGS  computer  code  HST3D  (Kipp,  1987) 
was  used  to  solve  the  equation  for  ground-water  flow 
and  solute  transport  using  backward-in-space  and 
backward-in-time  finite-difference  approximation 
methods.  Mass-fractional  concentrations  and  total 
fluid  pressures  were  simulated  using  the  model.  The 
reader  is  referred  to  the  original  report  that  documents 
the  computer  code  (Kipp,  1987)  for  a  complete  discus¬ 
sion  of  the  numerical  methods  and  model  code. 
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SPECIFIC  CONDUCTANCE,  IN  MICROSIEMENS  PER  CENTIMETER  AT  25  DEGREES  CELSIUS 
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Figure  10.  Changes  In  specific  conductance  during  the  9-day  aquifer  test  at  Walsingham  Park  Test 
Production  Well  and  Walsingham  Park  Monitor  Wells  1  and  2. 
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Conceptual  Model 

The  hydrogeologic  system  of  the  study  area  is 
conceptualized  as  a  multilayered  aquifer  system 
encompassing  various  permeable  zones  separated  hy 
confining  and  semiconfining  units  extending  from  the 
surficial  aquifer  system  downward  through  the  Suwan¬ 
nee  Limestone  of  the  Upper  Floridan  aquifer  (fig.  2). 
Ground-water  flow  is  considered  to  be  lateral  in  per¬ 
meable  zones  and  vertical  in  confining  and  semicon¬ 
fining  zones.  Ground-water  quality  changes  from  fresh 
to  saline  with  depth  as  ground  water  becomes  increas¬ 
ingly  dense  with  higher  concentrations  of  dissolved 
solids.  Several  assumptions  were  made  regarding  the 
conceptualized  hydrogeologic  system  to  effectuate 
numerical  simulation: 

1 .  The  aquifer  system  is  assumed  to  act  as  an  equivalent 

porous  medium, 

2.  All  aquifer  layers  are  homogeneous  and  isotropic, 

3.  Initial  conditions  are  hydrostatic, 

4.  Initial  fluid  density  is  constant  in  each  model  row, 

and 

5.  The  water-quality  profile  is  laterally  homogeneous 

for  the  entire  model  area. 

Design  of  Numerical  Model 

The  numerical  model  subdivides  the  conceptual¬ 
ized  hydrogeologic  system  from  the  ICU  through  Zone  B 
at  the  test  site  into  an  organized  network  of  grid  blocks, 
each  with  assigned  fluid  and  matrix  properties  repre¬ 
senting  a  500-ft  thick  section  underlying  the  test  site. 
The  model  was  constructed  by  using  a  cylindrical- 
coordinate  grid  of  124  columns  and  12  rows  (fig.  11), 
and  was  designed  to  simulate  the  responses  of  UZA 
and  LZA  to  brackish-water  development.  The  columns 
are  variable  in  width  with  radial  spacing  expanding 
logarithmically  from  0.3  ft  at  the  test  production  well 
up  to  a  maximum  width  of  328  ft  extending  out  to 

11.500  ft.  Radial  widths  extending  to  12,500  and 

14.500  ft  were  also  simulated,  but  the  additional  simu¬ 
lation  area  in  the  horizontal  direction  had  no  affect  on 
model  results.  The  12  model  rows  were  grouped 
together  as  five  distinct  hydrogeologic  layers.  Layer  1 
corresponds  to  model  rows  1-3,  is  175  ft  thick,  and  rep¬ 
resents  Zone  B  of  the  Upper  Floridan  aquifer.  Layer  2 
corresponds  to  model  rows  4  and  5,  is  50  ft  thick,  and 
represents  LZA  of  the  Upper  Floridan  aquifer.  Layer  3 
corresponds  to  model  row  6  and  represents  a  discrete, 
25 -ft  section  within  Zone  A  with  an  increased  occur¬ 
rence  of  green  clays.  This  section  was  modeled  to 


determine  the  hydraulic  properties  of  the  unit  separat¬ 
ing  UZA  and  LZA.  Layer  4  corresponds  to  model  rows 
7-9,  is  150  ft  thick,  and  represents  UZA  of  the  Upper 
Floridan  aquifer.  Layer  5  corresponds  to  model  rows 
10-12,  is  100  ft  thick,  and  represents  the  ICU. 

Initial  and  Boundary  Conditions 

Initial  conditions  refer  to  pressure  and  fluid  prop¬ 
erty  distributions  established  prior  to  the  initiation  of 
the  model  simulation.  Hydrostatic  conditions  are  estab¬ 
lished  in  HST3D  from  a  distribution  of  dissolved  solids 
and  fluid  densities  for  each  model  row.  Values  for  dis- 
solved-solids  concentrations  and  fluid  densities,  based 
on  water-quality  parameters  defined  from  samples 
taken  during  drilling,  were  used  to  establish  initial  con¬ 
ditions  (table  2).  The  model  calculates  the  starting  fluid 
pressures  based  on  the  user-specified  fluid  property  dis¬ 
tribution  for  all  model  rows  as  well  as  at  model  bound¬ 
aries.  Under  these  hydrostatic  conditions,  the  aquifer 
response  to  pumping  is  determined  by  model  simulated 
changes  in  pressure  relative  to  the  initial  hydrostatic 
conditions  (Anderson  and  Woessner,  1992). 

Boundary  conditions  delimit  the  model  domain 
by  constraining  the  lateral  and  vertical  extent  of  the 
simulated  flow  system  at  the  model  periphery.  The 
inner  cylindrical  boundary  represents  the  production 
well  and  was  defined  by  allocating  well  information  to 
the  various  rows  that  represent  the  open-hole  interval 
of  the  pumped  well.  Model  rows  not  open  to  the 
production  well  were  assigned  a  no-flow  condition, 
whereas  the  outer  cylindrical  boundary  was  defined  as 
a  specified-pressure  condition  (fig.  11).  The  assigned 
pressure  distribution  along  the  vertical  plane  of  the 
outer  boundary  was  calculated  by  using  water-level 
and  dissolved-solids  concentration  data  obtained 
during  drilling  of  the  test  production  well.  The  lower 
model  boundary  was  defined  as  a  specified-pressure 
boundary  equivalent  to  the  pressure  exerted  by  a 
column  of  stratified  saltwater  and  freshwater  in  the 
overlying  strata  determined  from  the  density  of  water 
samples  acquired  during  drilling  (fig.  11).  The  outer 
and  lower  specified-boundaries  were  used  to  simulate 
an  infinitely  extensive  aquifer  while  eliminating  the 
modeling  of  indefinite  regions  outside  the  study  area. 
The  upper  model  boundary  was  assigned  a  no-flow 
boundary  coinciding  with  the  top  of  the  ICU.  Because 
no  drawdown  occurred  in  either  the  SAS  or  the  ICU 
monitor  wells  during  the  aquifer  test,  a  no-flow  bound¬ 
ary  condition  at  the  top  of  the  ICU  was  deemed  appro¬ 
priate.  However,  row  1 1  (corresponding  to  the 
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Figure  11.  Model  grid,  simulated  open-hole  intervals,  and  boundary  conditions  used  in  the  HST3D  simuiation. 
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Table  2.  Values  for  model  parameters  used  to  establish  Initial  conditions  and  values  for  calibrated 
model  parameters 

[Ib/ft®,  pound  per  cubic  foot;  mg/L,  milligrams  per  liter;  ft^,  feet  squared;  ft/d,  feet  per  day] 


Initial  conditions  Calibrated  model 


Hydrogeologic  unit 

Density 

(Ib/ft^) 

Scaled-solute 
mass  fraction  / 
dissolved  solids 
(mg/L) 

Intrinsic 

permeability  (fF)  / 
hydraulic 
conductivity  (ft/d) 

Storage 

coefficient 

Effective 

porosity 

(percent) 

Intermediate  confining  unit 

62.38 

0.0089  /  250 

3.9x10-^^/0.0001 

3.1x10-"^ 

35 

Upper  Zone  A 

62.38 

0.0098  /  275 

OO 

' — 

b 

p 

i> 

2.6x10-"^ 

30 

Semiconfining  unit 

62.41 

0.0250  /  700 

b 

p 

— 

o 

4.3x10-^ 

35 

Lower  Zone 

62.61 

0.1786/5,000 

2.9x10-*^/ 750 

8.6x10-^ 

30 

(Row  6) 

Lower  Zone  A 

62.72 

0.2679  /  7,500 

2.9x10-^/750 

8.6x10-^ 

30 

(Row  5) 

Zone  B 

63.65 

1.0/28,000 

3.9x10-^2/i 

2.6x10-^ 

30 

®  Input  parameters  varied  between  the  two  rows  representing  Lower  Zone  A. 


open-hole  interval  of  WPMW-3)  was  used  as  an 
indictor  for  the  upper  limit  of  hydraulic  parameters  that 
would  simulate  a  no-response  condition  in  the  ICU. 
Although  the  lack  of  a  drawdown  in  the  ICU  made  it 
impossible  to  accurately  estimate  hydraulic  properties 
of  the  unit,  model  parameters  in  rows  representing  the 
ICU  were  adjusted  to  eliminate  drawdown  during  the 
simulation. 

Model  Parameters 

Model  parameters  including  hulk  matrix  proper¬ 
ties,  fluid  properties,  and  well  information  were 
defined  for  the  five  hydrogeologic  units  simulated 
(table  2).  Although  some  adjustments  were  made  to 
the  parameters  during  calibration,  initial  values  were 
carefully  selected  based  on  data  collected  during  this 
investigation  and  previous  studies  (Hickey,  1982, 

1989;  Knochenmus  and  Swenson,  1996;  Yobbi,  1996). 
Each  model  parameter  is  discussed  in  detail  below: 

I.  Bulk  matrix  properties. — Bulk  matrix  properties 
include  effective  porosity,  hydraulic  conductivity 
(as  a  function  of  intrinsic  permeability),  storage 
coefficient  (as  a  function  of  matrix  compressibil¬ 
ity),  and  dispersivity.  Bulk  matrix  properties 
(except  dispersivity)  for  each  hydrogeologic  layer 
are  listed  in  table  2.  Because  values  of  dispersivity 


are  poorly  known  within  the  study  area,  a  conser¬ 
vative  approach  was  used  by  setting  longitudinal 
and  transverse  dispersivity  to  zero. 

2.  Fluid  properties. — Fluid  properties  include  viscos¬ 

ity,  density,  compressibility,  temperature,  and  dis- 
solved-solids  concentration.  Dissolved-solids 
concentration  is  expressed  as  a  scaled-solute  mass 
fraction  ranging  from  0  to  1  in  HST3D.  A  value  of 
0  represents  pure  freshwater  whereas  a  value  of  1 
represents  water  with  the  highest  dissolved-solids 
concentration  (28,000  mg/L)  in  the  model.  Values 
for  density  and  dissolved  solids  for  each  hydro- 
geologic  layer  are  summarized  in  table  2.  Fluid 
viscosity  was  held  constant  for  the  entire  model  at 
1.0x10  centipoise.  Likewise,  fluid  temperature 
and  fluid  compressibility  were  held  constant  with 
values  of  77  °F,  and  3.03x10"^  in^/lb,  respectively. 

3.  Well  information. — Well  information  needed  for 

model  input  includes  the  withdrawal  rate  of  the 
pumping  well  and  is  distributed  equally  over  the 
three  rows  representing  the  pumped  zone.  The  test 
production  well  open  to  UZA  was  simulated  in 
model  rows  7-9.  The  simulated  pumpage  was  dis¬ 
cretized  with  a  constant  withdrawal  rate  assigned 
to  each  of  the  three  rows.  The  total  pumpage 
assigned  for  all  three  rows  for  the  aquifer  test  sim¬ 
ulation  equaled  220  gal/min.  Pumpage  was  varied 
thereafter  from  200  to  700  gal/min  during  the 
brackish- water  development  simulations. 
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Simulation  of  Aquifer  Test 

Numerical  analysis  of  the  9-day  aquifer  test  data 
provided  a  method  for  calculating  hydraulic  properties 
of  various  hydrogeologic  units  at  the  test  site.  The 
numerical  model  was  calibrated  hy  simulating  the 
aquifer  system’s  response  to  pumping  UZA  during  the 
aquifer  test.  The  model  analysis  of  the  9-day  aquifer 
test  consisted  of  meeting  two  calibration  objectives 
within  acceptable  limits:  (1)  simulate  the  drawdown 
and  recovery  of  water  levels  in  selected  wells  mea¬ 
sured  during  the  test,  and  (2)  simulate  the  overall 
response  of  the  system  to  pumping. 

The  first  objective  involved  matching  the  initial 
drawdown  part  of  the  aquifer  test.  The  initial  139  min¬ 
utes  of  the  aquifer  test  was  used  because  nearly  all 
drawdown  effects  occurred  during  this  time  and  moni¬ 
tor  well  drawdowns  could  be  separated  from  tidal 
effects.  Tidal  effect  corrections  were  made  by  calculat¬ 
ing  the  tidal  efficiency  for  the  monitor  zone  and 
adjusting  drawdown  and  recovery  data  accordingly. 
Tidal  efficiency  is  the  ratio  of  the  water-level  ampli¬ 
tude  measured  in  a  well  to  the  oceanic  tidal  amplitude 
(Ferris  and  others,  1962).  Tidal  data  were  collected  at 
a  gage  located  on  the  intracoastal  waterway  at  Indian 
Rocks  Beach,  2  mi  west  of  the  study  site  (fig.  1).  The 
calculated  tidal  efficiency  for  Zone  A  averaged  10  per¬ 
cent  for  ebbing  tides  and  16  percent  for  flooding  tides. 
After  tidal  corrections  were  made  to  the  observed  data, 
the  corrected  drawdowns  were  compared  to  simulated 
drawdowns. 

Model  parameters  were  adjusted  until  a  reason¬ 
able  match  was  derived  where  simulated  drawdowns 
closely  paralleled  the  observed  drawdowns.  A  reason¬ 
able  match  was  considered  less  than  0.10  ft  difference 
between  simulated  and  observed  drawdown  data.  Final 
calibrated  values  for  hydraulic  conductivity  were  1 8 
ft/d  for  UZA,  750  ft/d  for  LZA,  1  ft/d  for  Zone  B, 

IxlO  ’^  ft/d  for  the  ICU,  and  10  ft/d  for  the  semiconfin- 
ing  unit  separating  UZA  and  LZA.  Final  calibrated 
values  for  storage  coefficient  were  3.1x10'^  for  UZA, 
8.6x10'^  for  LZA,  2.6x10“^  for  Zone  B,  3.1x10"^  for 
the  ICU,  and  4.3x10’^  for  the  semiconfining  unit  sepa¬ 
rating  UZA  and  LZA.  Calibrated  model  parameters 
were  further  checked  by  comparing  the  observed 
recovery  data  of  the  aquifer  test  to  the  simulated 
recovery  water  levels  (fig.  12).  Values  of  the  corrected, 
observed  drawdown  and  recovery  and  simulated  values 
are  provided  in  table  3. 


Table  3.  Simulated  and  corrected  ebserved  drawdown  and 
recovery  for  selected  monitor  wells  during  the  9-day  aquifer  test 

[UZA,  upper  zone  A;  LZA,  lower  zone  A;  NA,  data  not  available; 
all  values  are  in  feet] 


Zone 

and 

weii 

Corrected 

observed 

drawdown® 

Simulated 

drawdown 

Corrected 

observed 

recovery'^ 

Simu¬ 

lated 

recovery 

UZA 

WPMW-1 

2.69 

2.66 

2.79 

2.66 

LZA 

WPMW-2 

0.29 

0.30 

0.36 

0.31 

UZA 

Botanical  Well 

NA 

0.59 

0.54 

0.60 

^Corrected  drawdown  for  first  139  minutes  of  pumping. 

'^Corrected  recovery  for  first  136  minutes  of  recovery  in 
WPMW-1  and  Botanical  Well,  and  first  192  minutes  of  recovery 
in  WPMW-2. 

The  second  objective  involved  matching  the 
long-term  pattern  of  drawdown  and  recovery  for  the 
length  of  the  9-day  aquifer  test.  As  figure  9  shows, 
after  the  initial  139  minutes  of  pumping,  observed 
drawdowns  stabilized  (excluding  tidal  effects)  for  the 
remainder  of  the  9-day  test  until  the  pump  was  turned 
off  and  recovery  effects  occurred.  A  reasonable  match 
was  accomplished  when  simulated  drawdowns  func¬ 
tioned  in  a  similar  fashion  with  drawdowns  decreasing 
after  the  initial  drop  in  water  levels  and  stabilizing 
throughout  the  simulation  until  the  recovery  part  of  the 
test  (fig.  13).  It  should  be  noted  that  the  simulated 
hydrographs  are  superimposed  over  the  uncorrected 
field  hydrographs.  No  attempt  was  made  to  simulate 
the  variations  in  water  levels  due  to  tidal  effects. 

Sensitivity  Analysis 

Sensitivity  analysis  was  used  to  determine  the 
effect  of  each  parameter  on  the  calibrated  model 
results.  Input  parameters  were  individually  varied  over 
a  reasonable  range  of  values  to  determine  the  sensitiv¬ 
ity  of  the  model  to  that  particular  parameter  over  the 
entire  simulation  period.  Ranges  were  based  on  values 
from  previous  investigations  discussed  in  the  hydro¬ 
geologic  framework  section  of  this  report.  Changes  in 
certain  parameters  had  no  effect  on  the  model, 
whereas  others  noticeably  affected  the  model  results. 
For  example,  a  one  order  of  magnitude  increase  and 
decrease  in  the  value  of  matrix  compressibility  in 
model  rows  representing  both  UZA  and  LZA  had  no 
effect  on  the  model  results.  Furthermore,  a  one  order 
of  magnitude  increase  and  decrease  in  the  intrinsic 
permeability  of  model  rows  representing  the  ICU  had 
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Figure  12.  Simulated  and  observed  water  levels  during  drawdown  and  recovery  periods  at  selected  wells 
during  the  9-day  aquifer  test. 
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no  effect  on  the  model  results.  Conversely,  the  largest 
changes  to  the  model  results  occurred  from  changes  in 
intrinsic  permeability  values  assigned  to  both  UZA  and 
LZA  (figs.  14a  thru  14d).  Model  rows  representing 
UZA  were  increased  and  decreased  by  one  order  of 
magnitude,  whereas  model  rows  representing  LZA 
were  changed  by  a  factor  of  2  and  0.5.  Additionally,  the 
value  of  matrix  compressibility  assigned  to  model  rows 
representing  the  ICU  were  increased  and  decreased  by 
one  order  of  magnitude  (figs.  14e  and  14f).  An  increase 
in  the  matrix  compressibility  of  the  ICU  resulted  in  a 
slight  response  in  only  UZA.  Furthermore,  this 
response  occurs  in  the  early  part  of  the  drawdown 
graph,  but  the  overall  magnitude  of  drawdown  in  UZA 
did  not  change.  The  model  was  generally  sensitive  to 
changes  of  intrinsic  permeability  and  insensitive  to 
changes  of  matrix  compressibility;  therefore,  more 
confidence  was  placed  on  estimates  of  transmissivity 
than  on  estimates  of  storativity. 


Model  Limitations 

Results  from  model  simulations  outlined  in  this 
report  are  approximations  of  the  actual  hydrogeologic 
system.  Confidence  in  model  results  is  dependent  on 
the  assumptions  outlined  earlier  in  this  report  and  the 
extent  that  the  hydrogeologic  system  is  simplified  for 
mathematical  representation.  Simplification  of  the 
conceptual  model  is  necessary  to  simulate  a  very  com¬ 
plex  natural  system;  however,  the  extent  to  which  the 
system  is  simplified  represents  a  source  of  error  in 
model  results.  The  major  limiting  factor  in  the  model 
simulations  is  that  the  hydrogeologic  system  was  sim¬ 
ulated  as  a  five-layered  system  with  isotropic,  homo¬ 
geneous,  and  porous-medium  individual  layers.  The 
actual  hydrogeologic  system  may,  however,  be  aniso¬ 
tropic,  heterogeneous,  and  possess  properties  of  sec¬ 
ondary  porosity  due  to  dissolution  and  fracturing. 
Although  alternate  combinations  of  model  parameters 
may  provide  similar  results  to  those  outlined  in  this 
report,  the  model  incorporated  the  best  estimates  of 
the  unknown  parameters,  the  local  geology,  and 
hydrologic  conditions. 
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EXPLANATION 

•  OBSERVED  WATER  LEVELS 
+  SIMULATED  WATER  LEVELS 


Figure  13.  Simulated  and  observed  water  levels  in  selected 
monitor  wells  during  the  9-day  aquifer  test. 
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WATER  LEVEL,  IN  FEET  ABOVE  SEA  LEVEL 


WPMW-1 


WPMW-2 


EXPLANATION 


.  OBSERVED  WATER  LEVELS  ~  MATRIX  COMPRESSIBILITY 

+  SIMULATED  WATER  LEVELS  k  INTRINSIC  PERMEABILITY 


Figure  14.  Sensitivity  of  simuiated  water  ieveis  to  changes  in  seiected  modei  input  parameters. 
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This  report  presents  the  results  of  a  study  to 
evaluate  the  hrackish-water  resources  in  the  Upper 
Floridan  aquifer  underlying  west-central  Pinellas 
County,  Florida.  A  production  test  well  and  four  moni¬ 
tor  wells  were  constructed  at  a  location  in  Walsingham 
Park  in  1996-97.  The  wells  were  sampled  for  water- 
quality  parameters  and  investigated  using  borehole 
geophysical  methods  and  hydraulic  testing.  The  data 
collected  were  analyzed  hy  using  numerical  tech¬ 
niques  to  delineate  the  hydrogeologic  framework  at 
the  test  site. 

Collected  data  indicate  that  a  multilayered  aqui¬ 
fer  system  consisting  of  two  permeable  zones  under¬ 
lies  Walsingham  Park.  One  zone  extends  from  about 
135  to  275  ft  below  land  surface  and  corresponds  to 
UZA  of  the  Upper  Floridan  aquifer.  Another  perme¬ 
able  zone  extends  from  about  300  to  350  ft  below  land 
surface  and  corresponds  to  LZA  of  the  Upper  Floridan 
aquifer.  Analysis  of  water-quality  data  indicates  that 
LZA  may  not  provide  water  of  consistent  quality 
because  of  rapidly  increasing  dissolved-solids  concen¬ 
trations  occurring  around  350  ft  below  land  surface. 
UZA  was  deemed  more  appropriate  for  production 
pumping  and  was  tested  during  a  9-day  aquifer  test. 

A  numerical  model  was  developed  to  simulate 
the  9-day  aquifer  test  performed  at  the  Walsingham 
Park  test  site.  The  model  was  calibrated  to  the  draw¬ 
down  and  recovery  water  levels  of  selected  monitor 
wells  at  the  test  site.  Final  calibrated  values  for 
hydraulic  conductivity  were  18  ft/d  for  UZA,  750  ft/d 
for  LZA,  1  ft/d  for  Zone  B,  1x10'^  ft/d  for  the  ICU, 
and  10  ft/d  for  the  semiconfining  unit  separating  UZA 
and  LZA.  Final  calibrated  values  for  storage  coeffi¬ 
cient  were  3.1x10“'^  for  UZA,  8.6x10“^  for  LZA, 
2.6x10"^  for  Zone  B,  3.1x10"“^  for  the  ICU,  and 
4.3x10"^  for  the  semiconfining  unif  separating  UZA 
and  LZA. 
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